Existence and topological stability of Fermi points in multilayered graphene 
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We study the existence and topological stability of Fermi points in a graphene layer and stacks 
with many layers. We show that the discrete symmetries (spacetime inversion) stabilize the Fermi 
points in monolayer, bilayer and multilayer graphene with orthorhombic stacking. The bands near 
k — and e = in multilayers with the Bernal stacking depend on the parity of the number of 
layers, and Fermi points are unstable when the number of layers is odd. The low energy changes 
in the electronic structure induced by commensurate perturbations which mix the two Dirac points 
are also investigated. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 



I. INTRODUCTION. 

The recent synthesis of monolayer graphite [l|, Q 
(graphene), the experimental ability to manipulate few 
layer samples 0, El 0| and the observations of quasi 
two dimensional behavior in graphite [6], have awaken 
an enormous interest in these systems. The conduction 
band of graphene is well described by a tight binding 
model which includes the it orbitals which are perpen- 
dicular to the plane at each C atom0, Q. This model 
describes a semimetal, with zero density of states at the 
Fermi energy, and where the Fermi surface is reduced 
to two inequivalent K-points located at the corners of 
the hexagonal Brillouin Zone. The low-energy excita- 
tions with momenta in the vicinity of any of the Fermi 
points have a linear dispersion and can be described by 
a continuum model which reduces to the Dirac equation 
in two dimensions J9. [lClf. what has been tested by re- 
cent experiments [J, LUJ . Fermi points have also been 
found in the modelling of the low energy band structure 
of multilayer systems both theoretically [H, EH and in 
experiments [lj, EH . A crucial issue for both theory and 
technology is the possibility of controlling the opening of 
a gap in the samples. From a theoretical point of view, 
the gap is related to the chiral symmetry breaking and 
mass generation, a classical -unresolved- problem that 
has been explored at length in planar QED [ll| QjJ • For 
the applications it is by now clear that opening a gap 
in monolayer graphene will be a difficult task and efforts 
are concentrated on multilayer structures [H, Ell EH . 

In this paper we analyze the stability of the Fermi 
points under small perturbations using very basic topo- 
logical concepts [2fJ. We find that the Fermi points are 
protected by the discrete symmetries (translational in- 
variance and space and time inversion) in the monolayer, 
bilayer AB and multilayers with ABC A ■ ■ ■ (rhombohc- 
dral) stacking. The stability of Fermi points in stacks 
with the Bernal stacking, ABAB ■ • • , depends on the par- 
ity of the number of layers. We also discuss the changes 



in the low energy and low momenta properties induced 
by commensurate perturbations which hybridize the two 
K points and partially break translation invariance. We 
will not analyze here in detail the effects of spin-orbit 
coupling, which my lead to additional changes at low 
temperatures [H Hi H3, H . 

The analysis reported here will be useful for the con- 
struction of continuum theories for long wavelength spa- 
tial perturbations, and for the study of the degeneracy 
and spectrum of the low index Landau levels in a mag- 
netic field. 



II. ELECTRONIC STRUCTURE AND 
STABILITY IN GRAPHENE. 

The Fermi surface (FS) is a central concept in con- 
densed matter that controls the low-energy physics of 
the systems. In a Landau Fermi liquid at T = 0, J. Lut- 
tinger [25] defined the FS of an interacting Fermi system 
in terms of the single-particle Green's function G(k,u>), 
as the solution of the equation G" 1 (fc, 0) = and showed 
that it encloses the same volume, equal to the fermion 
density n, as in the nonintcracting system. The robust- 
ness of the Fermi liquid idea has been understood re- 
cently in the context of the renormalization group where 
the Fermi and Luttinger liquids are seen as infrared fixed 
points[H, In recent works [2^. l29j Volovik has em- 
phasized the idea of the topological stability of the Fermi 
surface as the origin of the robustness of the Fermi liquid 
and has suggested a classification of general fermionic 
systems in universality classes dictated by momentum 
space topology. A more recent proposal relates the sta- 
bility of Fermi surfaces with K-theory, a tool used to 
classify D-brane charges in string theory [3(j. The idea 
behind the topological stability is to study the zeroes of 
the matrix Gq (fc, cj) (free inverse propagator) that can 
not be lifted by small perturbations. Here we will analyze 
the stability of the Fermi points of single and multilayer 
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FIG. 1: (Color online) Direct lattice and unit cell for mono- 
layer graphene. 
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electronic excitations behave like massless Dirac fermions 
with relativistic spectrum E = ±|fc|. 

Under a fc-independent, translationally invariant per- 
turbation 



H{K X + k) 



a z k* + a* 
k + a —a z 



(3) 



where a = a x + ia yi the spectrum becomes E = 
±\/al + \k + a\ 2 and a gap 2|a z |is generated. This is 
consistent with Horava's general arguments [3(j, which 
show that Fermi points for Dirac fermions in two di- 
mensions are unstable. However, this is not the end of 
the story, since the existence of discrete symmetries can 
sometimes stabilize the Fermi loci. In the case of the 
graphene, this role is played by time-reversal T : t — > —t 
and spatial inversion / : (x,y) — > (— x, — y). The real- 
ity of the 7T orbitals implies that time reversal merely 
reverses K 



(4) 



whereas the spatial inversion also exchanges the two 
types of atoms 



FIG. 2: (Color online) First Brillouin zone and Fermi points. 
The vectors K'±,K" (K' 2 ,K'i) are equivalent to Kx (K 2 ). 



I$ A (K) = $ B (-K) , I$ B (K) =$a(-K) (5) 

Invariance under these symmetries imposes the following 
constraints on the hamiltonian 



graphene, where the discrete symmetries of the system 
play a principal role. Although we will restrict ourselves 
to perturbations that can be studied within the context 
of a single-particle effective hamiltonian, the extension 
to self-energy induced perturbations is rather straight- 
forward and will be reported elsewhere [3l|. 

As shown in Fig. [TJ monolayer graphene consists of a 
planar honeycomb lattice of carbon atoms. Correspond- 
ing to the two atoms in the unit cell, one may define two 
Bloch wave functions to be used in a variational (tight- 
binding) computation of the spectrum 



t) 



A,B (1) 



where the sum runs over all the points in the direct lat- 
tice, i.e., t = niti +n2t2, {vat'Tb) are the positions of the 
atoms in the unit cell, and <&(r) is a real (7r-type) atomic 
orbital. As is well known, a simple tight-binding com- 
putation [7|, [8| yields a spectrum with two Fermi points 
located at K\ = — 2(fi/3 — 52/3 and K2 — —K\, where 
t.- L ■ gj — 2jrSij (see Fig. [5]). Near the two Fermi points, 
the hamiltonian can be linearized, and using appropriate 
units one finds 



H(Ki + k) 



k* 
k 



— kx@x 



k v a. 



(2) 



and H{—Ki + k) ~ —k x a x + k y a y . where k = k x + ik y 
and <7i are the Pauli matrices. Thus, the low energy 



T : H(K) = H*(-K) 
I : H(K) = a x H(-K)a x 



(6) 



Although these are useful properties that relate the 
hamiltonians at opposite values of K, what we need is 
a constraint on the form of H(K). This is obtained by 
combining time reversal with the spatial inversion 



TI : H(K) = a x H*(K)a x 



(7) 



implying Hu(K) = H 2 2(K). This enforces a z — in ([3|) 
and we see that no gap opens if the perturbation pre- 
serves the space-time inversion TI. 

This has an interesting topological interpretation, 
which extends the previous arguments to fc-dependent 
— but translationally invariant — perturbations. The low 
energy hamiltonian H{K\ + k) in ([2]) defines a map from 
the circle k x + k 2 = R 2 to the space of 2 x 2 hamiltonians 

H = h ■ a: 



Re 1 ' 



h y , h 



.) = i?(cos0,sin0,O) (8) 



Since Fermi points correspond to zeroes of the determi- 
nant —Det(H) = h 2 + h y + h 2 , a perturbation will be 
able to create a gap only if the loop represented by the 
map ([8]) is contractible in the space hamiltonians with 
non- vanishing determinants, which is just i? 3 — {0}. This 
is clearly the case, since iri(R 3 — {0}) = 7Ti(S' 2 ) = 0. On 
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the other hand, hamiltonians invariant under TI are rep- 
resented by points in i? 2 , and we have 



MR 2 ~ {0}) = ms 1 ) = z 



(9) 



This means that non-trivial maps such as the ones im- 
plied by ([2|) can only be extended to the interior of the 
circle by going through the origin, i.e., by having at least 
one zero. This precludes the creation of a gap. 

Note that the maps defined by the low energy hamil- 
tonian in the proximity of the two Fermi points 



H(±K 1 + ft) : k = Re w -*h x + ih y 

have opposite winding numbers N 
computed by the formula 



±Re 



±i.6 



(10) 



±1, which can be 



N : 



1 

4ni J 



■ / deTr{a z H- 1 d g H) (11) 
Jo 



The fact that the two Dirac points carry opposite 
charges suggests that they could annihilate mutually if 
brought together by a perturbation. Any external po- 
tential commensurate with the honeycomb lattice, which 
has a finite Fourier component at the wavevector G — 
K\ — K2, induces terms which hybridize the two Dirac 
points and it will lead to the possibility of a gap. We can 
compute all the possible perturbations which are compat- 
ible with the symmetries of the lattice. The most general 
(4 x 4) hamiltonian including perturbations at G = and 
G = Kx-K 2 = -Ki is: 



H = 



k* + Q\ 
fc + Qi 
Q2 Q4 

Qt Q3 





-ft* + Q\ 



<V3 

-k + Qi 




(12) 



where Qi = Qf + iQ\ transforms accordin g to the E2 
representation of the Cq v symmetry group [32| at the 
T(K = 0) point. Q 2 and Q3 belong to the E represen- 
tation of the C% v group at K± , and Q4 to the A\ rep- 
resentation of the same group (see [33| for notation). 
At this point it is worth noticing a point on notation. 
When grouping the hamiltonians attached to the two 
Dirac points (-Ki,a) into a 4-dimensional matrix it is a 
common practice to reverse the order of the sublattices 
(A,B) in one of the Fermi points in such a way that the 
4-dimensional wavcfunctions have the form 

If this is done the topological structure of the hamiltonian 
is messed up and the computation of charges becomes less 
clear. For this reason, we follow instead the convention 
in [2l|, where 

^ = {®Ki,A, ®Ki,B,®K 3 ,A, ®K 2 ,b)- 

This is also important if one tries to compare the analysis 
of the perturbations written in eq. (I12p with the ones 
produced by the different types of disorder [34], [H, [36| . 




FIG. 3: (Color online)Distortion caused by the condensation 
of A\ mode with real Q4 and new unit cell. 



The perturbation given by Qf and Q\ shifts the Dirac 
points, but does not open a gap. In fact, the only param- 
eter which opens a gap is Q4. When only Q4 is different 
from zero, the spectrum becomes E — ±y/\Q4\ 2 + \k\ 2 
and, for Q4 real, the deformation of the lattice is given 
in Fig. [3l where one can see that no point symmetry is 
broken. This shows, in particular, that invariance under 
space-time inversion is not enough to guarantee the sta- 
bility of Fermi points — translation invariace plays a cru- 
cial role: TI by itself makes de Fermi points individually 
stable, but they may still annihilate against each other 
in the presence of a perturbation that breaks translation 
invariance. A distortion of the type of Q4 can be induced 
by a substrate with a periodicity commensurate with the 
lattice, or by the effect of one layer on another when there 
is a lattice mismatch between them, as in samples grown 
on a substrate [!, ■ It can be responsible for the gap 
observed recently photoemission experiments 38]. It is 
worth noting that the perturbation denoted here Q4 has 
been studied in a graphene ribbon in [39( . 

When only Q2 or Q3 are different from zero, we find: 



£fe = ± 



\Q2,3 



± 



\Q 



2,3 



+ |ft| 2 . 



(13) 



The energy bands are represented in Fig. [4] for the par- 
ticular case Q 2 — 0, Q3 = 1. We can see that the spec- 
trum is the same as the one obtained in a simple model 
for a bilayer system pbT| . which will be discussed later. 
A complete analysis of the most general perturbation of 
the form (jT3J) will be given elsewhere (3lj . 

In the absence of time reversal symmetry, other per- 
turbations are possible, such as 



SH = Bxa z T z 



B 2 0-yTy 



(14) 



where a and r are Pauli matrices whose entries are the 
sublattice and K point indices respectively, and B\ and 
B2 transform like the z component of a magnetic field and 
are odd under time inversion. Note that the first term 
is the orbital part of the intrinsic spin orbit coupling in 
graphene [2l|, [22, 23, 24], and it opens a gap. The second 
term should appear in a general spin orbit hamiltonian 
which takes into account the coupling between the two 
K points. 
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FIG. 4: Energy bands along the line (0, k v ) for Q3 = 1, Q4 
0. 




FIG. 5: (Color online) Skecth of the three possible positions 
of a given layer with respect to the others in a graphene stack. 
Bernal stacking (1, 2, 1, 2, • ■ • , is described by two inequivalent 
planes, while orthorhombic stacking, 1, 2, 3, 1, 2, 3 ■ ■ ■ , requires 
the three inequivalent ones. 



III. ELECTRONIC STRUCTURE AND 
STABILITY IN MULTILAYERED GRAPHENE. 

The case of the multilayer is more interesting. We will 
concentrate on the ability of the TI invariance to pre- 
vent the creation of a gap. For the sake of definiteness, 
only staggered (ABA) and rhombohedral (ABC) stack- 
ing will be considered. The relative orientations of the 
ABC planes arc sketched in Fig. [5]- The two inequivalent 
atoms in layer n will be denoted (A n ,B n ). Our conven- 
tions are such that the couplings of an A n (B n ) atom to 
the three in-plane nearest neighbors are shaped as a Y 
(inverted Y), independently of n. Thus, the low energy 
limit of the "free" hamiltonian obtained by neglecting 
inter-layer couplings is block-diagonal, with 2x2 blocks 
given by @. 

The simplest model introduces interlayer hoppings t 
only between nearest neighbors. The resulting hamilto- 
nian for bilayer graphene in the vicinity of the K± Fermi 
point is 



H(k) 



( k* t 

k 

k* 

V t k Q 



(15) 



and the energy bands are given by (|13p with the replace- 
ment |<52,3 1 — * t. In the limit E <C t, one can ob- 
tain an effective hamiltonian [Z(| for the lowest energy 
bands. To this end, reorder the wavefunctions according 



(A 1 ,B 1 ,A 2 ,B 2 ) — > (A 2 ,Bi,Ai,B 2 ), so that in the new 
basis the hamiltonian becomes 



H(k) 



k* 
k 
k* t 
kOtO 



H\\ Hi 2 
H 2 i H 22 



(16) 



where is a 2 x 2 block. The identity 



Det(H - E) (17) 
= Det(H n - H 12 (H 22 - E)~ 1 H 2 i - E^j Det(H 22 - E) 

shows that, for E <C t, the substitution H 22 — E — > H 22 
reduces the computation of the lowest energy bands to 
the diagonalization of the 2x2 effective hamiltonian 



H eff = H n 



H\ 2 H 22 1 H 2 i 





k 2 



k* 2 




(18) 



This effective hamiltonian involves only the atoms 
(A 2 ,B\), which arc not linked by t and give rise 
to bands with zero energy at the Fermi points. 
Since (A 2 ,Bi) are interchanged under spatial inversion 
r ^> — r, the combined TJ-invariance imposes a con- 
straint H eff (k) = <T x H eff *(k)a x identical to ©. This 
implies (fc) = H 22 ^(k), which shows that no gap can 
open. According to (fTTj) the topological charge for the K\ 
Fermi point is +2 and, by time reversal invariance, the 
charge for —K\ is —2. Thus, as in the case of monolayer 
graphene, the Fermi points are stable under perturba- 
tions that preserve TI and translation invariance. For in- 
stance, a perturbation like trigonal warping [40] changes 
the off diagonal elements in eq. ljTSj) . k 2 — > k 2 + v 3 k* and 
splits the Fermi point of charge Q = +2 into three Dirac 
points away from the K point, and charge Q = +1, and 
another Dirac point at the K point and Q = — 1, but the 
total charge is conserved and no gap opens. However, a 
perturbation hybridizing K\ and — K\ or one breaking 
TI might lead to a gapped system with no Fermi points 
at all. A physical example is provided by the experiment 
described in [l8| where a gap is controlled by changing 
the carrier concentration in each layer. 

This analysis can be easily generalized to multilayer 
graphene with rhombohedral stacking. This type of stak- 
ing includes the links (B\ —A 2 , B 2 — A3, . . . , £?jv_i — An) 
and the effective hamiltonian, which involves only the un- 
linked atoms (Ai,Bn), is given by 



n 



1 



N-l 





k N 



k* N 





(19) 



The topological charge for the K\ (—K\) Fermi point is 
(-N). As the point group for multilayer graphene 
with rhombohedral stacking is D%d, which contains the 
inversion /, the system is invariant under TI, which inter- 
changes (Ai, Bn), and the whole argument goes through 
as before. Thus we conclude that the Fermi points for 
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FIG. 6: (Color online) Energy bands close to e = for a stack 
with 4 layers (left), 8 layers (center), and 9 layers (right). All 
stacks have the Bernal stacking. The Fermi velocity is vf — 
1, t = 0.1 and a layer dependent shift has been included: e n = 
{0.02,0.01} (left), € n = {0.02,0.01,0.005,0.002} (center), and 
e n = {0.02,0.01,0.005,0.002,0} (right). 



multilayer graphene with rhombohcdral stacking are sta- 
ble against perturbations that respect TI and translation 
invariance. 

The situation is very different for ABA stacking. An 
A^-layer graphene stack is invariant under the spatial in- 
version / only for even N, where the point group is D^d, 
while it is D%h for odd N. 

The 2N eigenstates in a stack with N layers at k = 
can be divided into two sets: N states at the orbitals 
connected by the interlayer hopping t, and N states in 
the other sublattice sites of each layer. If we only consider 
the hopping t, the first N states acquire a dispersion [4l| . 
lying in the range — 2t < e < 2t. The other TV states are 
degenerate with e = 0. A perturbation compatible with 
all the symmetries of the stack is a layer dependent shift 
of the onsite energies. This shift can be arbitrary, except 
for the twofold degeneracy related with the equivalence 
between layers which are symmetrically placed around 
the center, e n = tN-n+i- This is illustrated in Fig.[B]. 

The results in Fig. [5] show a gap at half filling in the 
stack with four layers, and overlapping bands at all en- 
ergies for the stack with eight and nine layers. Note that 
the gap in the stack with four layers does not require the 
existence of an external electric field, which will break the 
equivalence of the layers at opposite sides of the stack. In 
all stacks with an even number of layers, the two bands 
which start at the onsite energy of layers n and N — n + 1 
are degenerate at k — 0. The effective 2x2 hamilto- 
nian describing the bands near these degeneracy points 
have off diagonal elements with a non trivial phase as 
in eq. (|18p . Hence, the degeneracy has topological charge 
Q = 2, and it cannot removed by perturbations com- 
patible with the symmetries of the stack. For even n, 
an explicit computation shows that these two bands will 



have curvatures with the same sign near k = 0. Thus, the 
corresponding degeneracy points do not represent Fermi 
points. For odd numbered layers, the two bands disperse 
in opposite directions away from k — 0, and the degener- 
acy points become stable Fermi points at the appropriate 
doping. Note that the symmetries of the system allow 
for a direct trigonal-like coupling between the two lay- 
ers, which will split this Fermi point and give rise to four 
Fermi points showing linear dispersion, as in the bilayer. 
On physical grounds, this coupling will be negligible, un- 
less the two layers are contiguous. 

The low energy bands in a stack with an odd number 
of layers also contain doubly degenerate states at k = 0, 
associated to the equivalence between layers at opposite 
sides of the stack. But in this case the inversion / is not 
part of the symmetry group D^h of the stack, no invari- 
ance under TI can be imposed and, as a consequence, 
the first homotopy group ~k\ is trivial. This means that 
no conserved topological charge exists. Hence, a gap may 
open at k = when other perturbations consistent with 
the symmetries of the stack are included. Concretely, 
a direct coupling between orbitals in the same sublat- 
tice in layers separated by an odd number of other lay- 
ers will open a gap. Such a coupling, between layers 
which are second nearest neighbors, has been proposed 
in graphite [42l |. 



IV. CONCLUSIONS. 

We have presented a classification of the bands at low 
momenta and low energy of graphene layers and stacks 
with many layers. 

Each Fermi point in single layer graphene is stable 
against perturbations which preserve the discrete TI 
symmetry, and which do not mix the two Fermi points. 
A magnetic field, for instance, induces a gap in the spec- 
trum, see eq. (|14p . This term arises from the discreteness 
of the lattice, and it should be of higher order than the 
minimal coupling which leads to the formation of Lan- 
dau levels. Combining this and dimensional arguments, 
we expect it to be B\ oc vf/Ib x (a/^s), where a is the 
lattice constant, and lg = yj (e_B)/(c$ ) is the cyclotron 
radius. Thus, for B ~ 10T, we have B\ « O.lmeV. 

We have also classified the long wavelength perturba- 
tions commensurate with the graphene lattice, which can 
hybridize the two K points. Some of these perturbations 
open a gap in the spectrum, while others shift the posi- 
tion of the Dirac points. We expect that their strength 
will decay like a power law with the wavelength of the 
distortion. 

The low energy and low momentum spectrum of stacks 
with many graphene layers depend on the stacking order 
and the number of layers. For the most common case of 
the Bernal stacking, we find that layer dependent onsite 
energies lead to Fermi points with double degeneracy, 
topological charge Q = ±2, and a parabolic dispersion 
in k. This situation will be stable in stacks with a large 
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(even) number of layers. In stacks with an odd number 
of layers, there is no conserved topological charge and 
this degeneracy will be broken by additional interactions. 
Stacks with rhombohedral order lead to degenerate states 
with a large topological charge, Q = N, which will give 
rise to the formation of a cascade of Fermi points slightly 
away from k — 0, with lower topological charges. 

The most likely origin of the inequivalence between 
layers is the charge accumulation at the layers close to 
the surface[43]. An induced doping of 10 10 — 10 12 cm~ 2 
gives rise to shifts in the local potential of 0.01 — O.leV, 
so that the splittings associated to this effect can be eas- 
ily measurable. We find that a true gap opens, in the 
absence of an external field which breaks spatial inver- 



sion, only in a stack with four layers and Bernal stacking. 
Finally, stacking defects, which break the equivalence be- 
tween pairs of layers, will also break the degeneracy of 
all bands at k = 0. 
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